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Abstract 

We investigate the resonant process of nuclear excitation by electron capture, in which a con- 
tinuum electron is captured into a bound state of an ion with the simultaneous excitation of the 
nucleus. In order to derive the cross section a Feshbach projection operator formalism is intro- 
duced. Nuclear states and transitions are described by a nuclear collective model and making use 
of experimental data. Transition rates and total cross sections for NEEC followed by the radiative 
decay of the excited nucleus are calculated for various heavy ion collision systems. 
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I. INTRODUCTION 



In the resonant process of nuclear excitation by electron capture (NEEC) a free electron is 
captured into a bound atomic shell with the simultaneous excitation of the nucleus. It is the 
nuclear physics analogue of dielectronic recombination (DR), where a resonant excitation 
of a bound electron occurs. NEEC is the time- reversed process of internal conversion (IC). 
The excited nucleus can then decay radiatively or by internal conversion. In the later case, 
a resonant inelastic electron scattering on the nucleus occurs. 



The NEEC recombination mechanism has been presented^ 
Several studies have been made concerning NEEC in plasmas 



or the first time in Ref. 
[2} or in solid targets {3], 4. 
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In [3j, cross sections of the process are calculated through a scaling procedure applied to the 
results of DR, considering that the two processes differ only in their excitation mechanisms. 
Following that, in Ref. 4J, similar estimates of NEEC cross sections are obtained for the same 
nuclei by using experimental nuclear rather than atomic data. A more theoretical approach 
is provided in Ref. through an explicit treatment of the electron-nucleus interaction in the 



capture process, following the theory used for the calculation of IC coefficients from Ref. 
In Ref. [5|, non-relativistic calculations are presented for the case of NEEC into bare ions 
channeling through single crystals. The results are in disagreement with the previous ones 
from [3, (J As NEEC has not been observed experimentally yet, neither in plasmas nor in 
the case of ions in crystals or colliding with electron targets, the magnitude of its cross section 
is in doubt. Similar discrepancies exist also in the case of theoretical calculations for the 
similar process of nuclear excitation by electron transition (NEET). NEET is a fundamental 
but rare mode of decay of an excited atomic state in which the energy of atomic excitation 
is transferred to the nucleus. This corresponds to the time-reversed bound state internal 
conversion. Unlike NEEC, NEET has been observed experimentally [3], in the same year in 
which direct evidence of the bound internal conversion 8| has been reported, thus opening 
a new period in which experimental precision allows atomic shell models to have regard 
for the internal structure of the nucleus. The permanent development of the experimental 
techniques and the enhanced possibilities of preparing bare ions and electron targets make 
the experimental observation of NEEC a reasonable goal for the foreseeable future. Thus 
theoretical calculations for NEEC occurring in scattering measurements are particularly 
useful, especially in finding candidate isotopes and transitions suitable for experimental 
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observation. 

In both NEEC and NEET, that are at the border-line between atomic and nuclear physics, 
electronic orbital energy is converted directly into nuclear energy. They offer therefore the 
possibility to explore the spectral properties of heavy nuclei through atomic physics experi- 
ments. Experimental techniques developed for scattering studies of electron recombination 
with atomic ions, e.g. experiments with stored |9|, [lCj or trapped HQ ions, can be applied 
to gain information on the nuclear structure of several nuclides which is hardly accessible by 
nuclear scattering experiments. Especially, NEEC is expected to allow the determination 
of nuclear transition energies, the study of atomic vacancy effects on nuclear lifetime and 
population mechanisms of excited nuclear levels. 

In this work we calculate total cross sections for NEEC followed by the radiative nuclear 
decay for various transitions in heavy elements. Particular interest has been payed for 
collision systems where experimental requirements for the observation of NEEC are likely 
to be fulfilled. In order to derive the cross section formula for the process we extended 
the Feshbach projector formalism developed and used for the DR [13J to account for the 
interaction of electronic and nuclear degrees of freedom. The electric and magnetic electron- 
nucleus interactions are considered explicitly and the nucleus is described with the help 
of a nuclear collective model jl4|. The dynamics of electrons is governed by the Dirac 
equation as required in the case of high-Z elements. This formalism is presented in Section 
ITU The calculation of the NEEC rates for the electric and magnetic transitions as well as the 
nuclear model are described in Section IIII1 Numerical results of the calculation are given in 
Section ITVl We conclude with a short Summary (Section EJ). The derivation of the magnetic 
interaction operator related to Section [Til as we U as some larger formulas involved in Section 
IIHI are given in the Appendix. In this work atomic units have been used unless otherwise 
mentioned. 



II. THEORY OF NUCLEAR EXCITATION BY ELECTRON CAPTURE 

In this section we derive the total cross section formula for the NEEC process followed 
by the radiative decay of the excited nucleus. We consider that the electron is captured into 
the bound state in the Coulomb field of a bare nucleus and that only the nucleus is decaying 
by emitting a photon. However, the derivation of the cross section can be extended to the 
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many-electron case in a straightforward way and also an additional radiative decay of the 
electron can be treated by applying the formalism. 



A. Decomposition of the Fock space by means of projection operators 

The initial state j^) of the system consisting of the nucleus in its ground state, the 
free electron, and the vacuum state of the electromagnetic field can be written as a direct 
product of the nuclear, electronic, and photonic state vectors: 

= \N,pm s ,0) = \N) <g> \pm s ) <g> |0) . (1) 

Here, p is the asymptotic momentum of the electron, m s its spin projection, and \N) the 
nuclear ground state. The state \^/ d ) formed by the resonant capture has the form 

|* d > = \N*, n d K d m d , 0) = \N*) ® \n d K d m d ) <g> |0) , (2) 

with n d , K d , and m d being the principal quantum number, Dirac angular momentum, and 
magnetic quantum numbers of the bound one-electron state, respectively. The excited nu- 
clear state is denoted by \N*). The final state |^/) of the NEEC process contains a photon 
with the wave number k and the transversal polarization a = 1,2, and the nucleus which is 
again in its ground state \N): 

\\ff f ) = \N,n d K d m d ,ka) = \N) <g> \n d K d m d ) ® \ka) , (3) 

\ka) = 4J0) . (4) 

Here, the at is a photon creation operator. The corresponding conjugate annihilation 
operator is denoted by a^ a . 

To clearly separate these states in the perturbative expansion of the transition operator, 
we introduce operators projecting onto the individual subspaces. Characterizing the state 
of the electron in the positive part of the continuous spectrum by the energy e rather than 
the momentum of the free electron, we write the projector P belonging to the first type of 
subspace as 

P = y de ^ \ae)(ae\ . (5) 



4 



For brevity we introduce the multi-index a to stand for all discrete quantum numbers of the 
total system. The projection operator of the subspace spanned by states of the type (jSJ) is 
written as 

Q = 5>><?|, (6) 
q 

with the cumulative index q introduced again to summarize all discrete quantum numbers 
describing the bound electron and the excited nucleus. The subspace of the state vectors 
containing one transverse photon is associated with the projection operator 

1 kcr 

Assuming corrections due to two- or more-photon states [ijj and due to the presence of 
the negative electronic continuum to be negligible, we postulate the following completeness 
relation: 

P + Q + R=l, (8) 
where 1 is the unity operator of the Fock space. 

B. The total Hamiltonian of the system 

The total Hamiltonian operator for the system consisting of the nucleus, the electron, 
and the radiation field can be written as 

H = H n + H e + H r + H en + H er + H nr . (9) 

The Hamiltonian of the nucleus H n is written in terms of the nuclear collective model [ijj] 
by using creation and annihilation operators of the collective modes, B' x and Bx^, 

ff B = 5> A 2V^, (10) 

where Q\ are the phonon frequencies. The Dirac Hamiltonian of the free electron is given 
by 

H e = ca-p+ {(3- l)c 2 , (11) 
whilst the pure quantized radiation field is described by 

H r = J2 UJ k a h a kv ( 12 ) 

Act 



Interactions between the three subsystems are described by the three remaining Hamiltoni- 
ans in eq. (jUJ). We adopt the Coulomb gauge for the electron- nucleus interaction because it 
allows the separation of the dominant Coulomb attraction between the electronic and the 
nuclear degrees of freedom, 

H en = [ d 3 r n -^- . (13) 
J \r e -r n \ 

In eq. (ITS]) . p n (f n ) is the nuclear charge density and the integration is performed over the 
whole nuclear volume. The interaction of the electron with the transverse photon field 
quantized in unit volume is given by 



ka ' 

Here, a is the vector of the Dirac a matrices (a x , a y ,a z ), and e*£ CT is the polarization vector 
of the photons. Similarly, the interaction of the nucleus with the electromagnetic field is 
given by the Hamiltonian 

H nr = ~Y,J d * r ^^)\j% (^ka^Ha + h - C ) > ( 15 ) 
ka 

where j n {r n ) is the nuclear current. 

Using the projection operators we can separate the Hamiltonian as follows: 

H = H + V (16) 

with 

H = PHP + QHQ + RHR , (17) 

V = H -H = PHQ + QHP + PHR + RHP 

+ RHQ + QHR. (18) 

In this way the effect of the nuclear potential on bound and continuum electron states is 
included to all orders. The individual terms in the perturbation operator describe transitions 
between the different subspaces. For example, PHQ describes the transition of the bound 
electron to the continuum, i.e. IC, and QHP accounts for the time-reversed process of 
IC, namely, NEEC. PHR and RHP are the lowest-order operators for photoionization and 
radiative recombination, respectively. QHR and RHQ account for the radiative excitation 
of the nucleus and radiative decay of the nucleus or of the electron in a bound state. 



C. Perturbation expansion of the transition operator 



The transition operator is defined as 

T(z) = V + VG(z)V , (19) 

where G(z) is the Green operator of the system given by 

G(z) = (z - H)- 1 . (20) 

Here, z is a complex energy variable. The cross section for a process can be expressed by 
the transition operator as follows: 

-(E) = % lim \(* f \T(E + ie)\%)\ 2 p f , (21) 



(Kl k Fi e^0+ 

with the and as final and initial eigenstates of Ho, respectively. This cross section is 
differential with respect to the angle Qk of the photon emitted in the process. Fi denotes 
the flux of the incoming electrons, and pf the density of the final photonic states. 

We use the Lippmann-Schwinger equation to write the perturbation series for T(z) in 
powers of V with the Green function Gq(z) of the unperturbed Hamiltonian H : 

T(z) = V + VG (z)V + VG (z)VG (z)V + ... . (22) 

Since the initial state of the NEEC process is by definition an eigenstate of P, and the final 
state is an eigenstate of R, we consider the projection RTP of the transition operator: 

RTP = RVP + RVG VP + RVG VG VP 

+ RVG VG VG VP + ... (23) 

Here and in the following we omit the argument z. The first term in eq. fT23l) does not 
contribute to the NEEC process. After a further analysis of the second term RVGqVP = 
RH nr QGoQH en P and inserting the spectral resolution (jSJ) of Q in the second order in V we 
arrive to 

(jff f \RTP\ift ) ~ V (N,n d K d m dj ka\H nr \q)(q\H en \NJm s ^) ^ 

q Z E q 

The energy E® denotes the unperturbed eigenvalue of the state \q). If we continue analyzing 
the perturbation expansion ( 1231) . the term of third order in V can be written as 

RVG VG VP = RH er PG PH en QG QH en P 
-\-RH er PG()PH er RGoRH er P 
+RHQG QHRG RH er P . (25) 



The first two terms do not contribute to the cross section of the considered NEEC process. 
The last term is decomposed as 



RH QG qQ H RG qRH er P 



RH er QG qQ H er RG qRH er P 



4~ RH er QG oQ H nr RG qRH er P 



-\-RH nr QGoQH er RGoRH er P 



-\-RH nr QGoQH nr RGoRH er P . 



(26) 



Here, the first two terms are not considered, as they describe recombination by radiative 
decay of the electron. The process incorporated in the third term of ( 1261) is not possible. 
The remaining last term, QH nr RG$RH er P ', accounts for the capture of the free electron by 
exchanging a virtual transverse photon with the nucleus. As we show in Appendix [Aj it can 
be approximated by QH magn P, where 



is the magnetic interaction Hamiltonian. 

We continue the expansion (|23|) of the T operator and only consider the terms that 
contain QH en P as the first step and RH nr Q as the final step. The contribution of order V A 
can be decomposed as 




(27) 



RVGqVGqVGqVP = RH nr QG QH ne PG PH ne QGoQH ne P 



(28) 



+RH nr QG Q(H er + H nr )RG R(H er + H nr )QG QH ne P . 



We rewrite the first term as 



RH nr QGoQH ne PGoPH ne QGoQH ne P — 
Y J RHnrQG \q){q\H ne PG Q PH ne \q'){q'\G QH ne P, 



(29) 



and consider the diagonal matrix element 



(q\H ne PG PH ne \q) = / dej^ 



(q\H ne \ae)(ae\H ne \q) 
z-E° 



(30) 



with E° defined by Ho\aE°) = E°\aE°). Using the equality 




(31) 
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it can be further decomposed into 

j dc ^2 ^\ H ^\ a£ )^ a£ \ H ^\ ( l) _ A£ NP ^pIC ^ 

J ^ z E 2 

AE^ = V / cfe £ ^MMBM , (33) 
rf = iirJ2\(q\H n e\aE°)\ 2 . (34) 

a 

The notation AE^ P was introduced to denote the Coulomb nuclear polarization correction 
to the energy of the state q and r* c for its internal conversion width. V denotes the principal 
value of the integral. 

In a similar manner, the second term of ( 1281) can be analyzed. It can be separated into 
the following four parts: 

RH nr QG Q(H er + H nr )RG R(H er + H nr )QG QH ne P = (35) 
RHnTQGoQH er RGoRH eT QGoQH ne P 
+RH nr QGoQH nr RGoRH nr QGoQH ne P 
-\-RH nr QGoQH er RGoRH nr QG()QH ne P 
-\-RH nr QGoQH nr RGoRH er QG()QH ne P . 

The first term describes the emission and reabsorption of a photon by the electron recom- 
bined into the bound state. Its diagonal matrix element possesses a real and an imaginary 
part: 

(q\H er RG RH er \q) = A£f - V/ ad . (36) 

AEq E describes the one-loop self energy correction to the bound state energy of the electron. 
The imaginary part is the radiative decay rate of the electronic state and vanishes in the case 
of electron capture into the ground state of the ion. The second term in ( |35i) describes the 
emission and a subsequent reabsorption of a virtual photon by the nucleus, and its diagonal 
matrix element reads 

(q\H nr RG RH nr \q) = AE* SE - ^ rad . (37) 

Here, AE^ SE is the nuclear self energy correction to the energy of the ion, and r^> rad stands 
for the radiative decay width of the nucleus in the state q. The last two terms of fl35|) 
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incorporate corrections to the intermediate state energy due to the exchange of a virtual 
transverse photon between the electronic and the nuclear currents. These corrections are 
neglected in our treatment as they are expected to be far less than the overall accuracy of 
experimental nuclear excitation energies. 

Continuing the expansion (1231) of the T operator, the matrix element of the intermediate 
state Green operator in fourth order contains terms of the form 

( q \H lGo H t G H lGo HM) = £ m^EMMEj^m , (38 ) 

q> H 

where the label i stands for ne, er and nr. We adopt the so-called isolated resonances 
approximation by taking only the diagonal matrix elements into account, i.e. we set q = q' 
in (I3"51) . This approximation is valid if the distance between neighboring resonances is large 
with respect to their total natural widths, which is the case in all systems we studied. 
Higher-order terms can be summed then as a geometric progression 

1 ^ — > £.11 . . 

z~^E° ^ X = z - m 1 - x (39) 

1 k=0 1 

with the dimensionless quotient 

x = ((q\H ne PG PH ne \q) 
z-E° V 

+ (q\H er RG RH er \q) + (q\H nr RG RH nr \q)) , (40) 



resulting in 



1 



z-E\- A£ 9 NP - A£f - A£ g NSE + ^Tf + |r£' rad 



(41) 



Thus, the infinite perturbation expansion introduces energy corrections and widths into 
the energy denominator of the lowest order amplitude (|24|) . The final expression for the 
transition amplitude of NEEC into states denoted by d and followed by radiative nuclear 
decay is then 

(■* f \Kr{z)p\*i) = 

sr-^ (N, n d K d m d ,ka\H nr \^ d )(^ d \H en + H magn \N,pm s ,0) 

^ z-E°-AE d +ir d ■ [ ' 

Here we introduce the notation AE d = AE^ P + AE% E + AE^ SE for the energy correction 
and T d = T d c + r^' rad for the total natural width of the excited state \d) = \N*, n d K d m d , 0). 
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D. Differential and total cross sections for NEEC 



Equation (|2T|) gives the differential cross section in terms of the matrix element of the 
projected T-operator. Neglecting the interference of neighboring resonances, and taking into 
account only a single state d, with the corresponding magnetic substates, the cross section 
for a given reaction pathway i — > d — > / is 



d<Jl 7^ f {E) = ^ (NfljMr n d K d m d ,ka\H nr \NjI d M Id ,n d K d m d ,0) 



2 



\(N d I d M Id ,n d K d m d ,0\H en + H magn \NiIiM h ,pm s ,0)\ 2 
X V-Ej-MW + W (43) 

Ni as well as Nf represent the nucleus in the ground state, while N d stands for the inter- 
mediate excited nuclear state. The angular momentum I and its projection Mj are used to 
denote the nuclear states. One has to perform an average over the magnetic substates of 
the initial state of the system and a summation over the final states if these are not resolved 
in a NEEC experiment. The total cross section is calculated by integrating over the solid 
angle flk of the photon emission and averaging over the direction of the electron: 



*—.,(*) % E E 2,2,'+ 1, E 2I ' + ! E I J 



dVt p 



M If aMi d m d y ' M h m 3 " A/} 

dQk 



\(NfIfM Jf , n d K d m d , ka\H nr \N^I d M Id ,n d K d m d , 0)| s 



(e - E d f + ir 2 



4 1 d 

2 



x | (N*I d M Id , n d K d m d , 0\H en + HmagnWiliM^pm., 0) \ 2 p f . (44) 

We denote the corrected energy of the intermediate state by E d = E d + AE d . By intro- 
ducing the notation 

yi^d = 2 7T \^ \^ 

n 9/9.7". -I- "H 



Mj.m s Mj d rn d 



X 



2(2/, + 1) 

/ dn p \ (N* d I d M Id ,n M , 0\H en + H^mM^rn,, 0) f Pl (45) 



for the electron capture rate 

2tt 



At** 



x 



21,, + 1 ^ 

M lf a M Id 

dQ k \(NfI f M If ,n d K d m d , ka\H nr \N d I d M Id ,n d n d m d ,0)\ 2 p f (46) 
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for the radiative transition rate, and 



(E - E d f + \T d 



ME ~ E d ) = - ^ ; 2 (47) 



<Ti^ f (E) = — r - n L d {E - E d ) . (48) 



for the normalized Lorentz profile and taking into account the relation FiPi = p 2 /(2n) 3 , the 
cross section formula can be written in the condensed form 

2tt 2 A^Y iT A 

p 2 r 

Determining the total cross section of the studied process requires the calculation of the 
transition rates Y n and A r , and the initial and final state energies. In the actual calculations 
we neglect the additional corrections AE^ P and AE^ . The integration of the cross section 
over the continuum electron energy gives the resonance strength S d for a given recombined 
state d, 

f 2lT 2 Ad-^fyi^d 

S d = dE — r n L d (E - E d ) . (49) 

J P d 

The natural width T d of the nuclear excited state is of the order of 10~ 6 eV. In this interval 
the values of p 2 as well as of the NEEC rate Y n can be considered constant. As the Lorentz 
profile is normalized to unity, 

' dE L d {E - E d ) = 1 , (50) 
the resonance strength can be written as 

2-7T 2 Ad^fyi^d 

S d =^ Ar J n . (51) 

p 2 r d 

III. RATES FOR ELECTRIC AND MAGNETIC TRANSITIONS 

In order to calculate the NEEC rate we have considered the matrix element of the electric 
and magnetic interactions between the electron and the nucleus. We write the wave function 
of the system as the product wave function of the electronic and nuclear states, 

H fi = (N*I d M Id \(n d K d m d \H en + H magn \pm s )\N hM^) . (52) 

The initial state continuum electronic wave function is given through a partial wave expan- 
sion [161 ]. 

\pm s ) = Y ^ ( n e)C ^ ^ j- mi m s mj |e«m) , (53) 



rn, 
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where e is the energy of the continuum electron measured from the ionization threshold, 
e = \Jp 2 c 2 + c 4 — c 2 . The orbital angular momentum of the partial wave is denoted by I 
and the corresponding magnetic quantum number mi and the phases A K are chosen so that 
the continuum wave function fulfills the boundary conditions of an incoming plane wave and 
an outgoing spherical wave. The total angular momentum quantum number of the partial 
wave is j = \k\ — ~. The Yi mi (Q e ) denote the spherical harmonics and the C (I \ j; mi m s m) 
stand for the vector coupling coefficients. 

For describing the nucleus we have used a collective model [ljj] in which the excitations 
of the nucleus are assumed to be vibrations of the nuclear surface. The surface can be 
parametrized as 

oo A 

r(o, v,t) = R (i+Y;Y. • ( 54 ) 

A=0 n=-x 

Rq denotes the radius of a homogeneously charged sphere and the time-dependent amplitudes 
a\n serve as collective coordinates. Using this parametrization and requiring that the charge 
is homogeneously distributed, the nuclear charge density can be written as 

Pn(r,t)=po9(R(e,(p,t)-r) , (55) 

with the constant average density p = Jj^fjr- m the first order in the collective coordinates 
the density p n can be approximated as 

p n (r, t) = Po 6(R - r) + Po R 6(R -r)]T a* Xfl (t)Y Xfl (n) . (56) 

A/i 

The electron-nucleus interaction Hamiltonian in eq. ffl3|) can be written using the multipole 
expansion as 

Hen = J2 E ^TTl YLl{Qe) / rf3r ^TT^(^)Pn(r„) . (57) 
L=0M=-L J r > 

The radius r e (r n ) denotes the electronic (nuclear) radial coordinate and Q e (Q n ) stands for 
the corresponding solid angle. We can make the simplifying assumption that the electron 
does not enter the nucleus, therefore we assume that r e > r n . According to Ref. [ljj], this 
approximation should not affect the results for the considered transitions. The Hamiltonian 
can then be written 

00 L Att 1 

H en = Yl E 2L + l Y LM{^e)^lQLM , (58) 
L=0 M=-L e 
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with the help of the electric multipole moments 

Qlm = J d 3 r n r^Y LM (n n )p n (r n ) . (59) 

After performing some angular momentum algebra one finds for the matrix element of 
H en for a given partial wave component 

(N* d I d M Id | (K d m d | H en | Km) | N i I i M Ii ) = 



LM 



(2L + 1) 



xC(J< l d L- -M It M Id M) C{j jd L- -m m d M) C (j d L j; ± M . (60) 
RLj d j stands for the radial integral given by 

/CO 
drr- L+1 (f Kd (r)f £K (r) +^ d (r)^(r)) , (61) 

where g(r) and f(r) are the large and small radial components of the relativistic continuum 
electron wave function 

V £KTO (r) = , (62) 

and the g d (r) and f d {r) are the components of the bound Dirac wave functions 

= ( '^V*™ ) (63) 

with the spherical spinor functions Q™. For a given multipolarity L, the NEEC rate for an 
electric transition thus reads 

/ 1 1\2 



(64) 



where 



B T {EL, Ii - I d ) = in ±— i \(N*I d \\Q L \\NI l )\ 2 (65) 



is the reduced electric transition probability. 

If we consider the charge density of the nuclear collective model from eq. ( l56i) . the ma- 
trix element of the Hamiltonian H en can be conveniently written in terms of the reduced 
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transition probability B j without imposing any constraints on the electron motion. In this 
case the electric NEEC rate is given by 

Y " } = (2TTip fi ° 2(1+1) ' B T (EL ' h "* Id)(2jd + 1] 

x EKL/C^ij-ioi) 2 , (66) 
where the electronic radial integral is 



R lL = l^L-i / drrL+2 + g Kd (r)g £K (r)) + 



o Jo 



POO 

+ R L +2 drr- L+1 (f Kd (r)f £K (r)+g Kd (r)9(rU) . (67) 
JRo 

The magnetic Hamiltonian in eq. (|2"7|) can be written using the multipole expansion as 

H magn = -a-A = -- c J2 ■ Y L A l(tt e ) I d 3 r n ^j n (f n ) • F^*(O n ) . (68) 

lm J r > 

We use again the approximation that the electron does not enter the nucleus. Then the 
Hamiltonian reads 



2L + 1 V L 

LM 



H magn = 7^-7 V ^r~^M LM a ■ Y L ^(Q e ) , (69) 



where the magnetic multipole operator is given by [3] 



M LM = - % -\Jj^l J d 3 rr L Y£l(n n ) ■ j n (r n ) . (70) 
Here, the vector spherical harmonics are defined as [3] 

?LL(Pe) = J2J2 C ( L1L ' U( 1 M)Y Lv (n e )e q , (71) 



v q 

where q = 0, ±1 and the spherical vectors e q are 

1 



; V2 

e = e z , 
1 



fix ^^y ) j 



^ K e x -ie y ). (72) 
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Introducing the expression for the vector spherical harmonics in eq. ( |69|) we obtain 



LM 



X 



^F i ,(n e )|--ic(LlL; J /l -M)(a x + ^a y ) 



+ C(L li;i/0 - M)a z + ;-j=C(£ 1 L\ v ~ 1 - Af)^ - J , (73) 

where a x , a y and a z are the Cartesian a matrices. The matrix element of the magnetic 
Hamiltonian then yields 

(N*I d M Id \(n d K d m d \ - a ■ A^Km^NUM^ = 

4rri ^(-l) M 7^^ T (iV*/ d M / jM LM |iVJ J M /i ) 

LMu 

X 



( - 7^C(^ 1 L;v 1 - M)(n d «; d rnd|r- (i+1) y Ll/ (fi e )(Q; cc + ia y )|e«m) 
, C(L 1 L;u - M)(n d K d m d \r-( L+ VY Lu {tt e )a z \eKm) 
+^= C ( L 1 L ^ - 1 - M)(n d K (i m d |r- (i+1 ^ Ll/ (O e )(a :i: - ia y )|eKm)j 



+C(L 1 L; v - M)(n d ^m d |r-( i+1 )r L! ,(fi e )a z |£fi:m) 

(74) 

Introducing the notations „ and for the electronic matrix elements we get 

(N*I d M Id \(n d K d m d \ - a ■ A\eKm)\NIiM h ) = 



Airi ^(_i)^ + M J £ + * i)3 C(I d h L;M Id - M h M)(N*I d \\M L \\NI t ) 



LMu 

1 



X 



C(L1L;i/1 -M)T+ !U + C(LlL;vO — M) T div 



+ ^=C(LlL;u -1 -M)T-\ . (75) 

The explicit form of the electronic matrix elements can be found in Appendix [B] For a given 
multipolarity L, the rate for the nuclear excitation by electron capture for a pure magnetic 
transition is then 

{m) 16tt 3 (L + l) Pi ^ 
^ " L(2L + 1)3 5 T (ML ' /j /d) 2- 

X 



£(--^=C7(L1L;i/1 -M)4 + C(IU;,0 - M) 

1 



Mu 

+ -1 -M)T- ) . (76) 
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All the nuclear information is contained in the reduced magnetic transition probability 

B T (ML, Zj - 7 d ) = -^[(iV^llM^OII^)! 2 , (77) 

whose value can be taken from experimental data or from calculations involving different 
nuclear models. Given the different parity of the electric and magnetic multipole moments 
a transition of a given multipolarity L is either electric or magnetic. We consider only 
the cases of transitions with a certain value of L and we neglect the possible mixing ratios 
between electric and magnetic multipoles of different multipolarities. 



IV. NUMERICAL RESULTS 



We calculate total cross sections and resonance strengths for NEEC followed by the 
radiative decay of the excited nucleus for various collision systems, involving electric E2 and 
magnetic Ml multipole transitions. We consider the cases of nuclear isotopes with low-lying 
nuclear levels for which the NEEC process is more likely to be observed experimentally. 

For the case of electric multipole transitions we have considered the + — > 2 + E2 transi- 
tions of the | 6 U, | 8 U, If Cm, ^ 4 Yb, ™Et, JfGd, JfGd, If By and ^ 4 Dy even-even nuclei. 
The quadrupole excitations of even-even nuclei are well described by the collective model. 
For the calculation of the NEEC rate both formulas from eq. (1641) and eq. (1661) have been 



used for a comparison. A further E2 transition | — > ~ of the odd 6@ 3 Dy nucleus has been 
investigated. For this case we have calculated the NEEC rate assuming that the electron 
does not enter the nucleus. The reduced transition probability B ]_ (E2) for the even-even 



(if; 
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20j , and in the case of the 



nuclei as well as the energies of the nuclear levels were taken from 

l63 Dy isotope, from 21|. The nuclear radiative rate was calculated according to the formula 



where A = E,M stands for the type of transition. The two reduced transition probabilities 
for the emission, respectively the absorption of a gamma ray are related through the formula 

B I (XL, I d - If) = ^f^B T (XL, I f - I d ) . (79) 

The width of the excited nuclear state is then 

T d = Af~* f + Af c (80) 
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where Af c is the IC rate of the state d, related to the NEEC rate through the principle of 
detailed balance, 



For the NEEC rate we need to calculate numerically the radial integrals Rh,j d ,j that enter 
eq. (j64p and eq. fl66l) . We use relativistic Coulomb-Dirac wave functions for the continuum 



electron and wave functions calculated with the GRASP92 package [23j considering a homo- 
geneously charged nucleus for the bound electron. The value of RL,j d ,j is almost the same 
whether we use Coulomb-Dirac radial wave functions or we take into account the finite size 
of the nucleus. The finite size of the nucleus has a sensitive effect on the energy levels of the 
bound electron. The energy level of the bound electron is calculated with GRASP92 and it 
includes quantum electrodynamic corrections. The first term of the sum of radial integrals 
in eq. fl671) is about 3 orders of magnitude smaller than the second term, 

„L+2 



o Jo 



drr L+2 (f Kd (r)f £K (r)+g Kd (r)g £K (r)) 



POD 

«i? L + 2 / drr- L+l (f Kd (r)f £K (r)+g Kd (r)gUr)) ■ (82) 



Here, the nuclear radius Rq is calculated according to the formula [24] 

R = (1.0793A 1/3 + 0.73587)fm , (83) 
where A is the atomic mass number. If we make the approximation 

P CO 

4 2 L - R o +2 / drr ~ L+1 (/*(r)/«(r) + 9 Kd (r)9Ur)) (84) 
Jo 



the NEEC rate is exactly the one in (1641) . calculated with the assumption that the electron 
does not enter the nucleus. For the numerical cases of the even-even nuclei the difference 
between the rates calculated with eq. (JMj) and eq. fl66|) are from under 1% up to 6%. The 
difference is larger for the capture into the s orbitals and it increases with the value of the 
atomic number Z. For the capture of the continuum electron into the 2s orbital of gg 8 Cm, 
the value of the rate calculated using the non-penetrating approximation is 6% larger than 
the one calculated with the collective model, in which the restriction on the electron motion 
is avoided. 

For the cases of the U isotopes and for gg 8 Cm, the capture into the K shell is not pos- 
sible due to the low energy level of the first excited nuclear state. For these 3 isotopes, 
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TABLE I: Electric NEEC rates and resonance strengths for various heavy ion collision systems. 
E cxc i s the nuclear excitation energy, E c is the continuum electron energy and is the width of 
the excited nuclear state. 



Isotope 


^exc(keV) 


£ c (keV) 


Type 


Orbital 


Y n (l/s) 


r d (ev) 




S(barn-eV) 


66 u y 


73.392 


10.318 


E2 


ls l/2 


1 8fi 
l.oO • 


ins 

1U 


A Qfi 
4.oO • 


i n- 

1U 


-8 


O.OO ' 


in- 2 

1U 




78.591 


11.350 


E2 


1*1/2 


9 99 

z.zz ■ 


in 8 

1U 


K 7/1 
O. 1 4 • 


i n~ 

1U 


-8 


/i 7n 

4. 10 ■ 


1 n — 2 

1U 


™Yb 


76.471 


4.897 


E2 


1*1/2 


1 78 
1. / • 


in 8 

1U 


A &A 
4.o4 • 


i n- 

1U 


-8 


Q 9^ 

y.zo • 


m-2 

1U 


^Gd 


123.071 


64.005 


E2 


ls l/2 


K R7 
0.0 1 • 


in 8 

1U 


9 ^1 
Z.D1 • 


i n~ 

1U 


-7 


9 an 
z.yu • 


in- 2 

1U 


IfGd 


88.966 


74.742 


E2 


2*i/ 2 


0.04 • 


in 7 

1U 


1 91 
1.Z1 • 


i n~ 

1U 


-7 


7 n7 
/ .u t ■ 


1 n— 4 

1U 


g| 6 Gd 


88.966 


74.896 


E2 


2pi/ 2 


117 
L.LI ■ 


in 8 

1U 


1 Q9 
l.oZ ■ 


i n~ 

1U 


-7 


9 9fi 

z.zo ■ 


i n — 3 

1U 


64 ba 


88.966 


75.680 


E2 


2p 3 /2 


1 fin 


in 8 

1U 


1 97 
L.Zl ■ 


i n~ 

1U 


-7 


o.l t ■ 


i n — 3 

1U 




80.660 


65.432 


E2 


2*1/2 


9 fin 
Z.oU ■ 


1U 


y.oy • 


i n~ 

1U 


-8 


fi 9Q 
O.Zo • 


1 n— 4 

1U 




80.660 


66.594 


E2 


2pi/ 2 


1 fin 
l.OU • 


in 8 

1U 


111 
i.ii ■ 


i n~ 

1U 


-7 


9 QQ 

z.yy • 


i n — 3 

1U 


6 i 2 Dy 


80.660 


66.492 


E2 


2P3/2 


9 1 fi 
Z.IO ■ 


in 8 

1U 


i n/i 

1.U4 ■ 


i n _ 

1U 


-7 


A 9^ 
4. ZD • 


m-2 

1U 




73.440 


58.212 


E2 


2*1/2 


y.io • 


i n6 

1U 


i fifi 
1.00 • 


i n _ 

1U 


-7 




i n— 4 

1U 


efDy 


73.440 


58.374 


E2 


2pi/ 2 


fi QQ 

o.yo • 


1U 


1 Qfi 

i.yo • 


1 n~ 

1U 


-7 


O.D4 • 


1 n~ 4 

1U 


efDy 


73.440 


58.272 


E2 


2p 3 /2 


Q AA 
y.44 • 


in 7 

1U 


1 

l.oO • 


i n~ 

1U 


-7 


1 91 
1.Z1 • 


i n — 3 

1U 


236 tt 
92 u 


45.242 


11.113 


E2 


2*1/2 


1 1 fi 
1.10 • 


in 8 

1U 


1 SO 

i.oy • 


1 n- 

1U 


-8 


8 7Q 

o. /y • 


m-3 

1U 


236 tt 
92 u 


45.242 


11.038 


E2 


2P1/2 


3.16 • 


i n9 
10 


4.19 • 


10 


-7 


i n r 

1.05 • 


1 n— 2 
10 


236 tt 
92 u 


45.242 


15.601 


E2 


2p 3 /2 


3.22 • 


10 9 


2.16 • 


io- 


-7 


1.56 • 


KT 2 


238 tt 
92 u 


44.910 


10.782 


E2 


2*1/2 


1.17 ■ 


10 8 


1.90 • 


io- 


-8 


8.90 • 


10~ 3 


238 tt 
92 u 


44.910 


10.706 


E2 


2P1/2 


3.20 • 


10 9 


4.25 • 


io- 


-7 


1.06 • 


10- 2 


238 tt 
92 u 


44.910 


15.269 


E2 


2P3/2 


3.27' 


10 9 


2.19 • 


io- 


-7 


1.56' 


io~ 2 


IfCm 


43.380 


5.500 


E2 


2*1/2 


2.32 • 


10 8 


3.42 • 


io- 


-8 


1.79 • 


10- 2 


IfCm 


43.380 


5.398 


E2 


2pi/ 2 


5.61 • 


10 9 


7.42 • 


io- 


-7 


1.91 • 


io~ 2 


i 8 Cm 


43.380 


11.018 


E2 


2P3/2 


5.42 • 


10 9 


3.60 • 


io- 


-7 


2.20 • 


io- 2 



recombination into the L shell of initially He-like ions is the most probable one. We re- 
gard the capture of the electron into a closed shell configuration as a one-electron problem, 
without the participation of the K-shell electrons. We consider for the continuum electron 
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a total screening, while the bound electron is described by wave functions for an extended 
nucleus. The electron interaction is included in the bound radial wave functions calculated 
with GRASP92 and it influences the results through the value of the bound energy and 
through the shape of the electronic radial wave functions. The change of the shape of the 
radial wave functions has a larger numerical contribution to the value of the NEEC rate 
than the change of the energy due to the electron-electron interaction. 

For the other cases of even-even nuclei, capture into the K shell is possible. For the g| 6 Gd, 
gg 2 Dy, and gepDy isotopes we have also considered the capture into the He-like ions. The 
width of the nuclear state in eq. ()80p contains then an extra term which accounts for the 
possible IC of the K-shell electrons. The capture rate into the 2p orbitals is in general one 
order of magnitude larger than the one for the capture into the 2s orbital. The NEEC rates 
and resonance strengths for NEEC followed by the radiative decay of the nucleus for electric 
transitions are presented in Table [J . The values of the NEEC rates have been calculated 
using the non-penetrating approximation. 

For the magnetic multipole transitions we consider the Ml transitions of the odd isotopes 
If Bo, ^ 3 Yb, iMn, gFe, gK, JfGd, *f Gd, 7 f Re and £fRe. Numerical results for these 
ions are presented in Table [Til The electronic radial integrals are calculated numerically 
using the same type of wave functions for the bound and continuum electron as for the 



electric transitions. The reduced ma gne tic transition probability B { 



of the nuclear levels are taken from 
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26 



28 



29 



30 



31 



32 



Ml) and the energies 



331 ] . Capture into the K 



shell is possible for all the chosen ions, except for the IfGd isotope. Capture into the higher 
shells is less probable, and already for the capture into the 2s orbital of e7 7 Ho, the NEEC 
rate is one order of magnitude smaller. We present also results for captures into the He-like 
ions of 64 5 Gd and g| 7 Gd. The largest resonance strength is the one for the capture into the 
Is orbital of 7Q 3 Yb, namely, S = 8.57 • 10 _1 barn-eV. This value is small in comparison with 
the DR resonance strengths, which are in the order of magnitude of 10 3 barn- eV. 

In Figure [T] we present the cross section for the capture of the continuum electron into the 
2py 2 and 2p 3 / 2 orbitals of the two studied even-even uranium isotopes. The cross sections for 
the capture into the 2^3/2 orbitals are larger than the ones into the 2pi/2 orbitals. Although 
the cross section values are in the order of thousands of barns, the width of the Lorentzian 
is given by the width of the excited nuclear state, which is in the order of 10~ 7 eV. This 
validates the use of the isolated resonance approximation in Section III CI 
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TABLE II: Magnetic NEEC rates and resonance strengths for various heavy ion collision systems. 
E cxc i s the nuclear excitation energy, E c is the continuum electron energy and is the width of 
the excited nuclear state. 



Isotope 


^exc(keV) 


£ c (keV) 


Type 


Orbital 


Y n (l/s) 


r d (ev) 




S(barn-eV) 


165 tr 
67 n ° 


94.700 


29.563 


Ml 


ls l/2 


4.0U • 


1U 


7 11 
1 .66 ■ 


1U 


-6 


/I 


in- 1 

1U 


173 Y K 
70 


78.647 


7.073 


Ml 


l«l/2 


9 ^0 


in 9 


9 A 1 


1U 


-6 


8 ^7 


in- 1 

1U 


if Re 


125.358 


42.198 


Ml 


lSl/2 


q 9n 
y.zu • 


in 9 


1 ^1 
l.Ol • 


i n- 

1U 


-5 


7 


in- 1 

1U 


ifRe 


134.243 


51.083 


Ml 


ls l/2 


8 89 
O.oZ • 


in 9 


1 .0 ( • 


1U 


-5 


p. 

o.uo • 


1U 


|Mn 


125.949 


117.378 


Ml 


l«l/2 


1 11 
1.66 • 


in 7 


1 7^ 


1 n~ 

1U 


-6 


^ 09 


1 n— 4 

1U 


57 F 


14.412 


5.135 


Ml 


l«l/2 


2.24 • 


-i r*6 

10 


1.25 • 


10 


-9 


o o o 

8.83 • 


10 


SK 


29.829 


24.896 


Ml 


l«l/2 


6.03 • 


10 6 


8.85 • 


io- 


-8 


1.10 • 


10~ 3 


IfGd 


60.008 


45.784 


Ml 


2si/ 2 


9.38 • 


10 8 


8.21 • 


io- 


-7 


2.62 • 


10" 3 


6 fGd 


60.008 


45.938 


Ml 


2pi/ 2 


1.04 • 


10 7 


7.84- 


io- 


-7 


3.05 • 


10" 4 


IfGd 


60.008 


46.722 


Ml 


2P3/2 


3.75 • 


10 6 


7.80' 


io- 


-7 


1.08' 


10~ 4 


64 7 Gd 


54.533 


40.309 


Ml 


2si/ 2 


1.42 • 


10 8 


3.17' 


io- 


-7 


1.34 • 


10- 2 


157pj 
64 *jCL 


54.533 


40.463 


Ml 


2pi/ 2 


1.59 • 


10 7 


2.61 • 


io- 


-7 


1.82 • 


10 -3 


IfGd 


54.533 


41.247 


Ml 


2p 3 /2 


5.81 • 


10 6 


2.56 • 


10- 


-7 


6.67' 


10~ 4 



V. SUMMARY 

In this article we present a versatile formalism for describing complex processes actively 
involving atomic electrons and nuclei and derive a total cross section formula for the process 
of NEEC. The cross section can be traced back to the calculation of IC rates and radiative 
transition rates. It also requires the knowledge of nuclear excitation energies and the bound 
state energy of the electron after capture. 

We derive NEEC rates for the case of both electric and magnetic multipole transitions 
using relativistic electronic wave functions. For the bound electron we use radial wave 
functions which take into account the finite size of the nucleus. The nucleus is described 
using a nuclear collective model. The nuclear part of the matrix element is written with the 
help of the reduced nuclear transition probability whose value is taken from experiment. 

We calculate numerically the cross sections for NEEC followed by the radiative decay of 
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45 



30 



15 - 




11.038 15.601 
10.706 E (keV) 15.269 



FIG. 1: NEEC cross sections for Uranium isotopes as a function of the continuum electron energy 

the excited nuclear state for various nuclei. Particular interest is payed to the cases that are 
more likely to be observed experimentally. The obtained resonance strengths are typically 
4 orders of magnitude smaller than the DR ones, which is due to the small width of the 
excited nuclear states. 

If the initial and final states for NEEC are the same as the ones for the process of radiative 
recombination (RR), quantum interference between the two processes occurs. The cross 
section for RR is much larger than the one of NEEC, therefore, the inclusion of interference 
terms is expected to increase the cross section by several orders of magnitude, making the 
experimental observation of NEEC feasible. Calculations about the magnitude of this effect 
are in progress. 
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APPENDIX A: 



In this Appendix we show that the term QH nr RGoRH er P can be approximated by 
QH magn P, where we introduced the Hamiltonian 

1 r* f j3 jn( r n) 



Hmagn = OL \ d T n — — , (Al) 

c J \r-r n \ 

describing the magnetic interaction of the nuclear and electric currents due to the exchange 
of a transverse photon. The replacement is valid in the case when the frequency of the 
virtual photon is negligible, or, in other terms, when its wavelength is large compared to the 
typical linear size of the total system. 

In this derivation we use the second quantized forms 

ab u 

x / ^ M)S (W* + Mf.) (A2) 

and 

of the electromagnetic interaction Hamiltonians. The <p a form a complete set of one-electron 
states, and the cl (c&) are electronic creation (annihilation) operators. The ip s and ip t denote 
nuclear states and the B\ and B t are the mode operators of the nuclear collective model like 
in fflUj) . Here we only label these operators by one index for simplicity. 

Substituting these operators into the matrix element of QH nr RGoRH er P, we get 

£ HB^ = . EEE j A U.i (A4) 

r r r abst ao> J J ^ ' 

x (q\B\B t J d 3 rj n ■ e^ap^e^ \r) 

x(r|4c 6 J d 3 r e 0l(r e )a ■ ?^ai^e~ lkre b (r e )\ae) . 

Note that only the photon creation term of (IA20 and the photon annihilation term of (1A3j) 
contribute. Introducing the notation 



1 fc 



M e ab (k, a) = —^-j d 6 r e <jf a {r e )a ■ ^e^'fofc) , (A5) 
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M£(k', a') = -i- y 1 J d?rj n ■ e^'^ (A6) 

for the electronic and nuclear interaction matrix elements and taking into account k = k', 
a — a', (IA4I) can be condensed as 

(q\H nr \r) (r\H er \ae) 



E 



Z — E r 



EEE / ^^|r>(r|4c^|p> M&( ^y g) ■ (A7) 

Evaluating the above expression at z = e^ + e™, which is equal to the sum of the initial state 
electron and nuclear energies, only the r state for which E r = e e a + e" + ck holds has to be 
considered. Taking into account the property 

E<3-^><»-?i,)=3-*-£^ (A8) 

(7 

of the transversal polarization vectors that holds for any pair of vectors a and b, we arrive 
to 

E = E / rf 3, te|B t BiC t Cl|p) _^L . (A9) 

z — E r t-** J el — e e n — ck 

r abst ° b a 

M abst (k) denotes the product of electronic and nuclear matrix elements summed over the 
polarization directions 

M abst (k) = "£M e ab (k,a)M: t (k,a) = j d 3 r e <P:(r e ) 

x / d*r n e^-V [j n . a- {]n ' k) & ' k) ) b (f e ) . 



M abst (0) = -7^7 / d 3 r e <P* a (r e ) I d d r n j n ■ a<f> b {r e ) . (A10) 



In the long wavelength limit k — > 0, this further simplifies to 

1 1 

(2tt) 2 k 

Applying the identity 

— = V- - m5{w) , (All) 
w + ie w 

where V implies the principal value integration, and adopting the long wavelength approxi- 
mation, the real part of (IA9I) turns into 



lin 45Z / d 3 k{q\B\B t c\c b \ae)M abst (k) 



+ o2 

abst 



1 1 

+ 



el — e e a — ck + ie e e h — e e a — ck — ie 



• (A12) 



24 



It can be rewritten as 



^{q\BlB t c\c h \ae) j ' d 3 r e (fr* a (r e )V magn (r e ,r n ; e e b - e e a )(j) b {f e ) (A13) 

abst 

in terms of the effective magnetic potential 

V ma9n (n;st-sl) = -^-lim l',f'k l',f'r„ (]„■<< ( '■>" ' » ( " ' » 



2(27T) 2 ke^0j J ' l \ Ja k 2 

X e ^„-fe) / 1 _ + 



— e£ — c/c + £^ — — ck — 



1 4vr f 00 „ /■ rt _ (j n -k)(a-k) 

dim / dfc / rf 3 r n j n ■ a - v A 



(2 7 r) 2 |f n -f e |^o t / 7 n y n k 2 

x sin(A;|f n -r e |) ( — 1 ] . (AI4) 

W-^a-cfc + ^e eg-e«-cA;-2e/ 

In the long wavelength limit, the fc-dependent part vanishes, and V magn (f e ; el — el) turns 
out to be 

V magn (f e ; 0) = --a [ d 3 r nl ^- . (A15) 
c y \r n -r e \ 

This is equal to the magnetic Hamiltonian in (1A1I) Q.E.D. It can also be shown that 
the imaginary part associated with the Dirac delta term in (IA11I) vanishes if the frequency 
e e a — e\ of the exchanged photon goes to zero. 
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A further calculation of the three electronic matrix elements that enter, in eq. (1741) . the 
formula of the magnetic Hamiltonian, gives the following expressions for the first term T div : 



It. I' 



T + ^_ 2zf _ ird -i./( 2L+1 )( 2 ^ + 1 )( 2 ^ + 1 ) ( L l A l B > \ ( L I a h> 



Arc 



s< I i 1 • 1 1 A r, fi 1 ■ 11 

*C [ l A - j;m+ - - - m )C [l B , - j d ; m d - - - m d 



2" 7 2 2 7 V 2 " 2 2 



, v 3 ,,„„_! / (2L + l)(2^ + l)(2f B + l) L ^ Z B 

V/ I o j I ,> \ 

I I 1 \ / 1 I I \ f°° 

xC [l A > ^jd,m d ---m d jC [l B -j;m+- - - mj j drr- L+1 g Kd (r)f K (r) , (Bl) 
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for the second term v : 



I A In' 



T o i{ A(2L + 1)(2/ A + 1) (2/^ + 1) ( L l A l B ,\ ( 

dt ' u V 4tt o o J y 

xC (l A ^ j;m- ~ 171 ) C \ l B' \ 3d\m d - ~ m d J / drr~ L+1 f K J i )u : ( r \ 



v m—\ \—m d 



■H(-l) 



„ M+ i J (2L + l)(2Z A + l)(2l B , + l) U Ib- i /a i B ' 

4?r 'oOO/Vz/m + i-i-mrf 



/ 1 1 1 \ / 1 1 1 \ Z" 00 

xC g m + 2 _ 2 m ) ° V B ' 2 Jd ' md + 2 ~ 2 V / drr ~ L+1 f^ r )9n{r) 

+ . ( _ 1)md _i . /(2L + 1)(2Z A , + l)(2fe + 1) ( Ll A , l B \ ( L I a' l B 




/yz/-i-m d m + i 
1 1 

,,, , ,,.fiH — 

2 ■"*' " 2 2 7 V 2 2 2 



1 1 1 \ / 1 1 1 \ Z" 00 

( l A > - H\ m d + - - - m d ) C [l B - j;m+ - - - m ) J drr~ L+1 f K (r)g Kd (r) 



h/M ,„ l(l _i / (2L + 1)(2Z^ + 1)(2Z B + 1) ilx-ls £ ^ 

47r ' / \ i/ i - m d m - 



xC ^U' ^ jd; % - ~ m d j C ^l B lj;m-^m) / drr L+1 f K (r)g Kd (r) , (B2) 
and for the third term T^: 




) „ (1+ i i /(2L+1)(2Za + 1)(2Zb' + 1) ( L l A l B , \ [ L l A l B > 



L " »0 / \i/m-± -m d - \ 




> { ' ( ^ 2 J; m _ ^ m ) ^ B ' 2 Jd ' ™ d + 2~ ~ ^ md 



/■oo 

/ drr- L+1 f Kd (r)g K (r) 
Jo 



+2i{ _ ird+ hJ(2L + l)(2l A , + l)(2l B + l) I L l A , l B 



An 








/ 1 I I \ / I 1 1 \ f°° 

xC Na' - U\ m d + 7> ~ 2 md ) C Y B 2 T,m ~ 2 2 m ) J drr ^ L+1 9^ r )fn{r) . (B3) 



Here Z^ and Z B are the orbital quantum numbers for the upper and lower two component 
spinors of the initial continuum wave function. For the wave function of the final bound 
state the prime indices l A i and l B > are used. For a given value k, the following relations hold: 



Ia 



1*1-2' 



K if K > , 

|/c| — 1 if k < , 
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K — 1 if K, > , 



K 



if k < . 



(B4) 
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